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Abstract 

The solution to the Bahtsky-Kovchegov equation is found in the deep saturation do- 
main. The controversy between different approaches regarding the asymptotic behaviour 
of the scattering amphtude is solved. It is shown that the dipole amplitude behaves as 
1 — exp {—z + Inz) with z = ln(r^(5^) (r -size of the dipole, Qs is the saturation scale) in the 
deep saturation region. This solution is developed from the scaling solution to the homo- 
geneous Balitsky-Kovchegov equation. The dangers associated with making simplifications 
in the BFKL kernel, to investigate the asymptotic behaviour of the scattering amplitude, is 
pointed out . In particular, the fact that the Balitsky-Kovchegov equation belongs to the 
Fisher-Kolmogorov-Petrovsky-Piscounov -type of equation, needs further careful investiga- 
tion. 
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1 Introduction 

The main objective of this paper is to find the solution to the non-hnear Bahtsky-Kovchegov 
equation [1, 2] in the saturation domain. There exist two solutions for dipole scattering amplitude 
N which were found in Refs. [3, 4, 5]. They can be presented in the form: 

N{Y,r]b) = 1 - e-^(^) (1.1) 

where 

z = ln{Ql{Y;b)r^) = agCY + In (rVr^ (6)) (1.2) 
and Qs is a saturation scale [6, 7, 8, 9]. Constant C in Eq. (1.2) is defined as 

C = 

and X is the BFKL kernel [10] 

X(7) = 2z^(l) 

where = dlnr{'~f)/d'j and r(7) is the Euler gamma function. 

The value for the critical anomalous dimension is determined by the equation [6, 9, 11] 

X ilcr) 
1 Tcr 

In the first solution [3] 

m -- 

while the second one has the form [4, 5] 

0(^) = z 

These two solutions lead to quite different approaches to the saturation boundary A/" — > 1. 
However, since the equation is a non-linear one, we cannot claim that only Eq. (1.7) survives at 
high energies. 

Eq. (1.6) is well accepted by the experts, while Eq. (1.7) is still considered by the experts as 
shaky, mostly because its derivation has not reached a stage of transparency as the first solution. 

In this paper we present (i) the derivation of both solutions in the framework of the same 
method; (ii) a simple explanation of both solutions; and (iii) the general form of approaching the 
unitarity limit. 



X ilcr. 

1 - 7c 



(1.3) 



- "0(7) - "0(1 - 7) 



(1.4) 



dx (7c 

d ^cr 



2C 



\ii{z) 



;i.5) 



:i.6) 



;i.7) 
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2 Solution to the Balitsky-Kovchegov equation (general 
approach) 



2.1 Equation 

The Balitsky-Kovchegov equation[l, 2] which we solve in this section, has a form: 

dN{r,Y;b) Cpag f £r' 



t/ ^"/^ ^ 2 1^/2 



(2.1) 



(^2N (^r',Y;b+ -{f-f')j - N {r,Y;b) - N (r' ,Y;b - -{f - f')^ N - f\Y;b - -f' 

where N {r,Y; b) is the scattering amplitude of interaction for the dipole with the size r and 
rapidity Y = ln(l/a;) {x is the Bjorken variable), at impact parameter b. 

It is useful to consider the non-linear equation in a mixed representation, fixing the impact 
parameter b, and introducing the transverse momenta as conjugate variable to the dipole sizes. 
The relations between these two representations are given by the following equations 

roo 

N{r,y;b) = rM kdk Jo{kr) N{k,y;b) ; (2.2) 
Jo 

N{k,y;b) = / - Mkr) N{r,y;b) ; (2.3) 
Jo r 

In this representation the non-linear equation reduces to the form [6, 13, 4] 

dN{k,y;b) ^ as{xm))my-b) - N'{k,y-b)) (2.4) 



dy 

and xiliO) operator defined as 



7(0 = 1 + ^ (2-5) 



where = ln(A;^i?^), and k is the conjugate variable to the colour dipole size and R is the 
size of the target. In this definition of the variable ^ we implicitly assume that b <^ R, and 
the amplitude N does not depend on b. The alternative approach for large values of the impact 
parameter was developed in Ref. [4], where the definition of variable ^ is quite different. 



2.2 Goal and assumptions 

Our main goal is to find the solution in the saturation region where r^Q^ ^ 1- Due to s- 
channel unitarity constraint [14] ^ the scattering amplitude in space representation (A^) should 

^See also Ref. [15] for an application of the Froissart boundary for hard processes 
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be less than 1. In the momentum representation this hmit means that N ^ We can obtain 
this estimate just using Eq. (2.3), and noticing that for r > 1/k the integral over r is small due 
to oscillating behaviour of Jq in Eq. (2.3). Therefore, for small k Eq. (2.3) reduces to 

N{k,y;b) = / -N{r,y;b) = - / dHl/r') N{r,y;b) (2.6) 

Assuming that N —1 for all r > 1/Qs (ln(l/r) < ln((5s)) we obtain from Eq. (2.6) that 

N{k,y;b) ^ Un(Ql/k') (2.7) 

Therefore we know the asymptotic behaviour of the amplitude, and we need to find how our 
function approaches Eq. (2.7). 

In the saturation region we expect the so called geometrical scaling behaviour of the scattering 
amplitude which was proven in Ref. [9], for such equations (sec also Ref. [13] for more general 
arguments and a more rigorous proof, and Ref. [12] for an observation that geometrical scaling 
behaviour could be correct, even in the part of perturbative QCD kinematic region). It means 
that N{k, y; b) is a function of the single variable 



ln(g^(z/,6)A') = as^^^ (Y -Yo) - ^ - m ; (2-8) 



where ^cr is a solution of Eq. (1.5) [6, 11, 13] and Yq is the initial rapidity. Function P{b) depends 
on impact parameter, but we will not discuss it here. 

Introducing a function 0(2;) we are looking for the solution of the equation in the form 

N(z) = ^J'dz'(l - e-*(^')) ; (2.9) 

Eq. (2.9) includes the geometrical scahng behaviour and leads to the asymptotic behaviour of 
Eq. (2.7). 

Our assumption that function is a smooth function, such that (p^z -C (pz <i>z where we denote 
02 = d(f)/dz and (pzz — d^(f)/{dzY is essential.. This property allows us to rewrite 

Jn 

{dzY ^ ^ ' ^ ' 

Eq. (2.10) means that we can use the semi-classical approach for the solution to Eq. (2.1) [4]. 

2.3 Reduction of the Balitsky-Kovchegov equation to the equation in 
one variable in the saturation domain 



Substituting in Eq. (2.4) N in the form of Eq. (2.9), and replacing F by i we obtain 

Xilcr) dN{. 
1 - Jcr dz 



«5 ^ = as ( x(l - /) m - N\z) ) (2.11) 



where / denotes / = d/dz = —d/d^ in Eq. (2.11). 

Differentiating both part of Eq. (2.11)with respect to z we reduce Eq. (2.11) to the form 



^e-*(^)0;(i) = fx{l-f)N{z)-N{z){l 



= ^(/X(l-/) - l)^(^) + N(z)e-^(^^ (2.12) 

An important property of function /x(l — /) — 1, is the fact that at small / it has an 
expansion that starts ^ from Since in our case / is operator / = ^, it means that the 
operator /x(l ~ /) — 1 contains the third and higher derivatives with respect to z. Therefore, 
one can see that the first term on r.h.s. of Eq. (2.12) is proportional to e~'^^^^ ( see Eq. (2.10)). 
Cancehng e"'^ on both sides of Eq. (2.12), and once more taking the derivative with respect to z 
we reduce Eq. (2.12) to the form: 



1 — 7cr {dzy ^ ' d(f)z (dz) 

Function dL/dcp'^ decreases at large values of the argument but has double pole singularities 
in all integer points (0^ = 1, 2, 3, . . . ). 

2.4 Two solutions. 

The existence of two solutions with sufficiently different forms of the dipole amplitude approaching 
its asymptotic value {N — 1 in space representation and N — |ln(Q^/A;^) in the momentum 
representation) can be seen directly from Eq. (2.13). Indeed, wc expect that (f){z) is large at large 
values of z and , therefore, we can neglect the term e^'^^^^ in Eq. (2.13). The first solution can be 
obtained from Eq. (2.13) assuming that dL/d(j)'^ gives a small contribution while (j)'{z) is large. If 
it is so Eq. (2.13) reduces to the simple equation 



Xilcr) 



1 - 'Jcr (dz) 

which leads to 



(2.15) 



7cr Z'^ 



^It should be stressed that the simphfied model for x(l — /) — f(^i_f^ does not have this property. This is an 
explanation why in Ref. [3] where this model was used, the solution was missed as we will discuss below. 
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at large z. The clearest case when we, indeed have this solution, is the simplification of the BFKL 
kernel which was considered in Ref. [3]. Namely, the kernel was taken as 



1 f 2 n,i 

asN, 



TT 



forr^g^^i. 



1-7 



(2.17) 
for r'^Ql > 1 ; 



instead of the full BFKL kernel uj{-f) = xil)- 

In this model function L is equal to zero and Eq. (2.16) gives the only solution in the saturation 
domain. 

For the full BFKL kernel there exists another possibility to have a solution: the term with 
dL/d(j)'^ compensates 1 in the r.h.s. of Eq. (2.13), and the l.h.s. is still small. 

Function being small for < 1 has a singularity at — > 1 namely 

— ; o for cpz 1 

d^z (1 - (t>zf 

Therefore, the first requirement leads to the equation 



(2.18) 



(1 - 0.)' {dzf 
For large z Eq. (2.18) has a solution 

<p{z) = z - \nz (2.19) 

which can be verified by explicit calculations. It should be stressed that (f){z) of Eq. (2.19) 
satisfies all conditions of a smooth function that has been used for the derivation of Eq. (2.13). 
Eq. (2.10) holds since (p^z ~ 1/^^ < 0^ 1. We can also check that the l.h.s. of Eq. (2.13) 
is proportional to 1/z^ and it can be neglected. 

We check the vahdity of the solution of Eq. (2.19) in more direct weay searching for the 
correction to Eq. (2.19) due to a violation of Eq. (2.10). The need to do this, arises from the 
appearance of a contribution of the order of 1/2;^ in 0^ which we cannot guarantee. Searching for 
such corrections we replace Eq. (2.10) by a new equation 



e 



= {{-4>zT - n{-(t>zr-^(t>zz) e-^(^' (2.20) 



{dzy 

Using Eq. (2.20) we obtain 

Xilcr) d'' (j) ^ (, _ dL{(f)z) d^(j) d^L{(j)z) d^ (j) ^ 

1 - ^cr {dzy \ J d(t)z {dzf {d(t)zY ^{dzf' ^' ' 
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instead of Eq. (2.13). Eq. (2.18) has the following form 

\ l± - 1 + I .2 22) 

(1 - <t>.f {dzf - '+(1 - M'\dzf ^'-''^ 

where we denote by 0o the solution of Eq. (2.19). Substituting 0^ = 0o,z + in Eq. (2.22) we 
obtain that 

dA(l)z 2 .c?Vov2 .2 



dz (1 - 00,.) ^dzf"^ ^ P ^^'^^^ 

One can see that Eq. (2.23) leads to A^^ oc 1/z^ and the solution is 

(i){z) ^ z - Inz + I ^ (po{z) + 0{l/z) (2.24) 

z 

Hence, Eq. (2.19) is a solution. 

It should be stressed that Eq. (2.19) is a solution to the homogeneous equation ( Eq. (2.13) 
with the l.h.s. equal to zero). Therefore, the dependence on Y in this solution stems only from 
the matching of this sohition to the linear equation at 5 < 0. As has been shown (see [6, 11, 12]) 
the solution of the linear equation behaves as 

N{z) = No exp((l-7e,)i) (2.25) 

We will discuss this matching in more details later. 

2.5 The complete solution and matching with the pQCD domain 

Assuming that in Eq. (2.13), 0^(i) is a function of 0(5) we can rewrite this equation in the form 

f xilcr) ^ dL^\ ^ ^ ^ _ ^_,(,) ^2.26) 



\ 1 — 7cr d(f)'^ J ^ dcf) 

Integrating Eq. (2.26) with respect to we reduce this equation to the form which gives the 
implicit solution for i as a function of 0. 



1 - 7cr 2 



= + e"^ - 00 - (2.27) 

where 0^o 0o are the initial conditions, namely, the value of function and 0^ at i = 0. 
There are a number of relations between them. The first one comes from the matching of the 
logarithmic derivatives at 5 = 0: 

1(1 -e-M 

= 1 - 7c. (2.28) 

iVo 
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The second relation originates from the matching of the logarithmic derivatives for dN/dz 
using the fact that the function N has the following form near to the saturation line:iV(^) — 
Noe^^~'^'''^'> ^ (see Eq. (2.25)), and that the geometrical scaling behaviour works even for negative 
z [12]. This relation is 

= (1 - 7cr) (l-e-^o) . (2.29) 

Unfortunately, with the initial conditions given by Eq. (2.28) and Eq. (2.29), we could not 
solve Eq. (2.27) analytically. Fig. 2.5 presents the numerical solution of this equation for positive 
5 as a function of the initial condition, namely, A^o (see Eq. (2.28) and Eq. (2.29)). 



100 




^ 80 
N 

9- 60 

II 

40 

9- 

20 



L , , , , , ^ 

20 40 60 80 100 120 

Z 

Figure 1: The numerical solution to Eq. (2.26) as a function of the initial condition Nq. Curves 
1, 2 and 3 correspond to Nq — 0.1,0.5,0.65, respectively. The lines are the asymptotic solution 
(f){z) — nz — \n.z -\-\n.\n.z with n = 1, 2, 3 for curves 1,2, and 3 , respectively. 

One can sec that the type of asymptotic behaviour depends on the value of A'o at ^ = 0. It 
has a simple explanation since dL{(f);^)d(f)z {n — If 0z at 2; = is smaller than unity we 

have the solution given by Eq. (2.19). However, if 2 > > 1 at 2; = 0, we have (t)z{z) — > 2 at 
large z. For 3 > 0^ > 2 the asymptotic behaviour is 3^; — In^;. 

It should be stressed that we cannot trust Eq. (2.27) for 2; — > since we cannot justify 
Eq.(2.10) in this region. In vicinity of 2; = this equation holds for the function N{z) rather 
than for the function (j){z). Nevertheless if A^o 1 0o is small as well and our approach could 
be justified. From Fig. 2.5 one sees that at small Nq we have the asymptotic behaviour given be 
Eq. (2.19). 
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3 Linearized equation in the saturation domain 



The key problem is the fact that the master equation is a non-hnear one. Therefore, generally 
speaking, we cannot conclude that the asymptotic behaviour of the solution can be written in the 
form (in space representation) 

N{z) = 1 - e-'f" - e-^ (3.1) 
where 0i is the solution of Eq. (2.16), while 02 is given by Eq. (2.19). 

However, the form of Eq. (2.9) shows us that the corrections are small in the region of large z 
and , we can therefore try to find the linear equation in this region for the function e""^. Denoting 
by S'(r, Y] h) we obtain the following equation for S from Eq. (2.1): 

^^'C-^"'") = (3.2) 



dY 

^ /p^^(s('-T;S) - s{r',Y;i-\(r-?))s{r-?,Y;i-\f 

The first solution comes from the region of integration r' ~ r 3> l/Qs for the second term 
on l.h.s. of the equation. In this region S <^ 1, and therefore, the non-linear term in Eq. (3.2) 
can be neglected. The linear equation has a very simple form 

dY ^ ~~ir k/Qi y-r'^\r'^ S[T,Y-h) = -as ln(Q, r )S[T,Y-h) (3.3) 
Substituting a new variable z (see Eq. (1.2)) Eq. (3.3) has the form: 

X(l-7cr) dS{z) 



1 - 7cr dz 

and solution to Eq. (3.4) is very simple, namely. 



-zS{z) (3.4) 



'''' 

The simple derivation of this solution makes it transparent. 

The second solution comes from quite a different region in integration in Eq. (3.2), namely, r' 
OT r — r' are much smaller than r and, basically, of the order of 1/Qs since our solution depends 
on z. Let us assume for the sake of presentation that 1/Qs ~ |r — r^l ^ r. For such small 
distances we can replace S{f—f', Y, b) by 1 {S{f—f', Y, b) = 1). Indeed, in pQCD region 1 — S* is 
the scattering amplitude and this amplitude is small. Therefore, the linear equation that governs 
the asymptotic behaviour of S has the following form: 

dS{r,Y) Cpas f d\' r 



dY TT 



? J (73^^ {Sir,Y) - 25(r-',y)|,_,,|^,/Q^<<,) (3.6) 
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In Eq. (3.6) we neglected the impact paramer dependence of S. 

One can recognize in Eq. (3.6) the BFKL equation, but with an overall sign minus in front of 
the r.h.s., and with the restriction that the second term is valid only if |r — f'\^ ^/Qs ^• 
The factor 2 in front of the second term on the r.h.s. of Eq. (3.6), arise from the fact that we 
have the same contribution from the kinematic region where r' ~ 1/Qs <S f. We can rewrite 
Eq. (3.6) in the form which is even closer to the BFKL equation by replacing the second term on 
r.h.s. of Eq. (3.6) by 

/ (7^^^^(-^^)|.-.'|.VQ.<<. - / (7^^ (^(-"'^) - ^^(-'^)) (3-7) 

In Eq. (3.7) we subtracted the region of integration l/Qs <^ r' <^ r or/and l/Qs <^ r' ^ r. 
Indeed, rewriting the kernel in Eq. (3.7) for r' < r in the form (after integrating over azimuthal 
angle) 

one can see that the first term has been taken into account in Eq. (3.6) since the region r' r gives 
the dominant contribution to this term. In this statement we assume implicitly that S{r\ Y) is 
steeply decreasing in the saturation region. We have to subtract the second term. We can replace 
r' by r in 5" since this function is large only at r' = r. Returning to the full kernel we obtain 
Eq. (3.8). 

Taking Eq. (3.8) into account we can rewrite Eq. (3.6) in the form 

dS{r,Y) Cpas f cPr' 



dV 



^ / 773W7^ (^(-'^) + {^(-'^) - ^S{r-,Y)}s.KL) (3.9) 



The second term on the r.h.s. of Eq. (3.9) is the kernel of the BFKL equation. 

Eq. (3.9) is very similar to the BFKL equation, but it has a negative sign in front and also an 
additional term, proportional to S (r, Y). 

Using the form 

S{r',Y) = 5(^) e-^^ (3.10) 

we can see that Eq. (3.9) can be reduced to Eq. (2.13). However, we do not need to take into 
account all the terms in Eq. (2.13). The main contribution comes from the first term on the r.h.s. 
of Eq. (3.9) and from the second term in the region of integration r' > r. Indeed, substituting 
Eq. (3.10) in the second term on the r.h.s. of Eq. (3.9) the relevant contribution appears as 



r 

Jr2 



dr 



12 




1 = (3.11) 



[ ^ {t^'^ - i) = m - ^{1 + 



1 - 0i 
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Using the variable z instead of Y in Eq. (3.9), we obtain the simple equation which contains 
both discussed solutions: 

X(l-7cr) dS{z) , I 1 (o^o\ 

Using S {z) — e~'^(^) we can easily reduce Eq. (3.12) to an algebraic equation for 0'^;, namely, 

^J,(.) = - + (3.13) 

with the solution 

= ^ + C ± ^{z + CY + 4C(l-z)) (3.14) 
where C is given by Eq. (1.3). At large z 

4>'}^^ - ^ (3.15) 

<t>'t^ - 1 - - (3.16) 
z 

The second term on the r.h.s. of Eq. (3.12) is very small if we substitute the solution (f)'}^'^ of 
Eq. (3.15) into Eq. (3.12). Therefore, this branch leads to 

S,{z) ^ e-^^ = e-r ^ exp [-^ (3.17) 

The branch ^'^.^"^ leads a small l.h.s. of Eq. (3.12) and therefore, the solution leads to as a 
cancellation of the first and second terms on the r.h.s. of Eq. (3.12). The solution takes the form 

S2{z) = e-'^' ^^-r dz'<t>'J-\z') ^ zexp{-z) (3.18) 

The solution S2{z) is the solution to the homogeneous equation ( see Eq. (3.9) with dS{r, Y) / BY = 
). The entire dependence of this solution on energy (Y) stems from the matching with the 
solution of the linear equation for negative z. 



4 General solution for approaching the unitarity bound- 
ary 

Here, we are going to find a general solution to Eq. (3.9). Wc can solve this equation which is 
an equation of the BFKL -type, using Mellin transform which we use for solution of the BFKL 
equation: 

S{z) = ^ / dfs{f)ef^ (4.1) 



11 



The main observation is that the function (r'^/r^)-^ is an eigenfunction of the BFKL kernel [10]. 
As we have seen (see Eq. (3.11)) the integration in the region r, > r leads to the eigenvalue 

/ {t~^ - l) = V^(l)-^(l-/) (4.2) 

For r' < r we have the following integral after integration over azimuthal angle : 



-z + ^ {tf-'- - 1) = -z + m)-m) (4.3) 

Collecting Eq. (4.2) and Eq. (4.3) we reduce Eq. (3.9) to the form: 

Cfsif) = ^ + Xif)s{f) (4.4) 

To obtain Eq. (4.4) we used the result that a multiplication by z translates into operator 
—d/df for the Mellin transform. 

Eq. (4.4) has the following solution 

s{f) = Ae^'^f^ (4.5) 

where 

X{f) = / d f {C f - xif)) = C^-2^(l)/ - lnr(l-/) + lnr(/) (4.6) 
Jo 2 

where constant C and x{f) ^-^e defined in Eq. (1.3) and Eq. (1.4) , respectively, while '0(1) is the 
Euler constant, which is equal to 0.577216. 

The general solution to Eq. (3.9), which satisfies the initial condition S{z = 0) = Sq, can be 
written in the form 

S{z) = So -L-/^cxp(X(/) + /;^) (4.7) 



a— too 



2m X{f) 



where a in the contour of integration is situated to the right of all singularities of the integrand 
[10]. 

At z = Eq. (4.7) leads to 

1 r dX 

S{z = 0) = So-- / — = So (4.8) 
12 



One can sec two important properties of the integrand in Eq. (4.7): (i) there are no singularities 
in it for / > 0, since the singularities of r(l— /) in the exponent cancel the singularities of X'[f) in 
the numerator of Eq. (4.7); and (ii) contribution at the pole / — > vanishes since the numerator 
is equal to zero in this point. This fact arises due to the cancellation of single pole and double 
pole contributions. This observation means that we correctly evaluated the integral of Eq. (4.3). 
There is a problem in that we can justify the second term in the r.h.s. of Eq. (3.9) only for 
\r — r'\ 1/Qg <^ r. However, at / — >■ the region of |r — r*| ~ r also contributes 
since ln(r'^/r^) ^ 1//. Strictly speaking we needed to subtract term 1// from Eq. (4.3). Our 
observation means that it was correct to keep this term in the evaluation of the large z behaviour 
ofS{z). 

The integral in Eq. (4.7) has two sources for the asymptotic behaviour: the saddle point at 
large / and the pole (singularities) contributions at / — > —n where n — 0,1,2 

To calculate the saddle point contribution at large values of / we replace all Gamma functions 
in Eq. (4.7) by the asymptotic expression. Therefore, the expression for X{f) has the form 



(4.9) 



X{f) ^ C ^ + 2/ ln(/) at / » 1 

The equation for the saddle point is 

Cf^-'J + 21n(/(+)) + ^ = (4.10) 

and 

fsV - ~ - '2Hz) (4.11) 
Therefore, the saddle point contribution has the form 



S^{z) = '^o^^exp(-lc(^ + 2 1n(^))2) (4.12) 



which is the same as our solution of Eq. (3.17), but obtained with better accuracy. 

Near the singularity / = — n the integrand is 



Closing the contour about the point j = —n, we see that the contribution is equal to 

poo rco 1 



We used here the integral representation of the Euler gamma-function (see Ref. [19] 8.314) 
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One can see from Eq. (4.14), that only the first singularity is essential since all others are 

suppressed at large values of z. 

However, they give more than the saddle point, and because of this, we have to sum all these 
contributions if we want to take into account the saddle point contribution as well. The sum gives 

S-^{z) = 1^(1 - exp(-e-)) (4.15) 
Finally the behaviour of the solution in the deep saturation region can be written as sum 

SasymAz) = {z) + S+ {z) (4.16) 

At large z Eq. (4.16) leads to 

{z) ^ Sorrr^e-' (4.17) 

This solution in the form of Eq. (4.16) has a very clear physical meaning. Each term is 
the suppression of gluon emission. Indeed, approaching the black disc limit, we expect that an 
emitted gluon could not survive, except at the impact parameters close to the edge of the hadron 
disc [16, 17, 18]. Inside of the disc, the dipole could elastically rcscattcr but cannot emit gluons. 
The term in Eq. (4.16) describes the probability of dipole interactions without emission of 
a gluon in the leading twist. All other terms correspond to such a probability but for higher 
twists contributions in the BFKL equation, since the poles at j — —n at n > 1 correspond to the 
contribution of higher twists to the BFKL Pomeron ^ (see for example Ref . [20] ) . 



5 Numerical check of the new solution 



Despite the analytical calculations that led us to the solution in the form of Eq. (4.17) it is 
necessary to check that the homogeneous equation ( our master equation with zero l.h.s. ) has 
a solution. The main reason for such a check, is that the analytical consideration is correct only 
for large ^ ^ 1, and the matching between negative z and large but positive z is still out of 
theoretical control. The second motivation for searching for a numerical solution is the result of 
the numerical simulation by Salam [24] who confirmed the solution of Ref. [3] rather than the 
solution of Eq. (4.17). 

We search for the solution of the following equation 



dN{y,q) 

dy 
as ■ Nc 

TT 



(5.1) 




{q - q') 



N{y,q') - 



q'^ + {q-q'f 



N{y,q) 



- N\y,q) 



^One of us (E.L.) thanks J. Bartels for very useful discussions on higher twist contributions to the BFKL 
equation, that he had with him a number of years ago. 
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This equation is the BaHtsky-Kovchegov equation but in the momentum representation, where q 
is a conjugated variable to r while Q is a conjugated variable to 6. For simplicity in Eq. (5.1) we 
put Q = 0. After integration over angles Eq. (5.1) reduces to the form 



dN{y,q) 

dy 
as ■ Nc 

TT 



(5.2) 



coo 



jV(j;,9") 

\q' ^ — g^l 



N{y, q^) 

ql 2 |g/ 2 _ g2| 



+ s_ 



N{y,q^ 



q'2 ^4g/4 + q4 



N\y,q') 



We wish to solve this equation assuming the geometric scaling behaviour of the solution, namely, 
this solution is a function of the only one variable 



z = \n(Ql{y,b)/e) = as^^^ (Y - Y^) - ^ - /3(6) ; 

J- Tcr 

We can express Eq. (5.2) in terms of single variable z using the following expressions: 



(5.3) 



1. q 
2. 



and dq"^ 



— e ^ dz 



d_ 

dy 



l-7cr OZ 



Thus, finally we get following non-linear equation: 

dN{z) 



dz 

1 - 7cr 
Xilcr) 



dz' 



N{z') 



N{z) 



+ 



N{z) ■ e 



y'(2e-^')2 + (e--)2_ 



- N\z) IhA) 



As we have mentioned Eq. (2.9), we can find the scattering amplitude N{z) using function 
(l){z). We are going to demonstrate numerically, that (f){z) — — {z — log 2;) (Eq. (3.18)) minimizes 
the right hand side of Eq. (5.4), i.e. it can be considered as solution for homogeneous equation: 







7c 



Xilcr) 



dz' 



N{z') 



N{z) 



— e~ 



— e~ 



+ 



N{z) 



^/(2" 



-z'\ 2 



+ e- 



-z\2 



- N\z)\{^.^) 



In order to perform required calculations we have to expand Eq. (3.18) to the negative values 
of z. It is well known [11, 12, 13] that at the negative z, the scattering amplitude behaves as 
N{z) = A^oe^^"^''''-'^. We also require continuation of N{z) and N'{z) at 2; = 0. 

We also require the matching of the two solutions at ^ = namely, the value of the amplitude 
N{z) and its derivative N'{z) should be equal at 2; = 0. 



N{z) 



ATo + \ - /?-exp 



z < ; 



- {z' - log [z' + i-2^oa->.) i 



] )]) dz' ^ > ; 



(5.6) 
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Parameter /?, actually, accumulates all information which is beyond the precision of our solu- 
tion. As we can see from Fig. 2-a variation of this parameter changes the value of the r.h.s. of 
the equation (see Fig. 2) that is calculated using Eq. (5.6). One can see that the r.h.s. of the 
equation vanishes if we choose a value of (3. It means that (p = z — Inz really is a solution to the 
homogeneous equation which fulfills the correct boundary condition. 




Fig. 2-a Fig. 2-b 

Figure 2: In this plot we demonstrate results of numerical calculations of evaluation Eq. (5.5) for 
different values of (3 with initial condition: N{z = 0) = 0.2 (Fig. 2-a) and N{z = 0) = 0.4 
(Fzg. 2-b). 

We hope that this numerical calculation will dissipate the lingering doubts that the homoge- 
neous equation has a solution, which approaches the asymptotic solution for high energy ampli- 
tudes. 



6 Beyond the Balitsky-Kovchegov equation 

In Ref. [12] a new approach beyond the mean field approximation is proposed, based on the 
statistical interpretation of the high energy behaviour of the scattering amplitude. The main 
result of this paper can be understood in the following way. 

The scattering amplitude for dipole-target interaction can be viewed as the average of the 
solution to the Balitsky-Kovchegov equation e with the saturation scale Qs, over the saturation 
scale. The average saturation scale is given by the equation [21] 

dlnQl _ _ Xilcr as (1 - 7cr) X'^lcr) 

dY """l-Tcr 2 ln\l/al) ^ ^ 

and averaging should be performed with a Gaussian weight of variance 

= casY/\n\l/as) (6.2) 
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Therefore the dipole amphtude is equal to [21] 

A{Y,r) = I dpN{z,p)e-^ (6.3) 

where z = Ql) and p = IniQ^Ql). 

This averaging drastically changes the asymptotic behaviour of the first solution N oc e ac . 
It transforms into 

N oc e~2c A{z,Y) ex e ^ (6.4) 
Using Eq. (6.2) one can see that Eq. (6.4) leads to the following asymptotic behaviour at fixed 

A{z,Y) (X exp(^-|^CMn3(l/a/) (6.5) 

On the contrary averaging does not dramatically change Eq. (4.17) leading mostly to redefi- 
nition of the saturation scale, namely, the new scaling variable is equal to 

C ^ z - ^a' (6.6) 

or, in other words, instead of saturation scale Qs given by Eq. (6.1) we have to introduce a 
saturation scale 

g; = a;e--^«exp(4c^r) (6.7) 

In new variable the asymptotic behaviour of the scattering amplitude looks as 

Therefore, this solution shows the geometrical scaling behaviour while the first one displays 
no such scaling behaviour. 




7 Summary 

We found that the Balitsky-Kovchegov homogeneous equation (see Eq. (2.1) with the l.h.s. is 
equal to zero) has a solution. This solution is determined by the singularities of the BFKL kernel 
in the anomalous dimension (Mellin conjugated variable to In(r^) where r is the dipole size). This 
solution has the following form 

N{r^) = l-Cexp{-z + \nz) (7.1) 

where z = ln(r^/i?^) where R is the arbitrary scale. C is an arbitrary constant that we cannot 
calculate. Setting = V/Q'g{Y) ( Qs is the saturation scale) on physical grounds we obtain that 
the solution to the Balitsky-Kovchegov equation has the form: 

N{r^) = 1- Const ex.p(^-\n{r^Ql) + \n\n{r^Ql)) (7.2) 
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where Const cannot be determined within our accuracy. 

Since it depends on the singularities of the BFKL kernel it is very dangerous to make any 
simplification of this kernel. Indeed, in two approximate expressions for the BFKL kernel that 
have been suggested (see Refs. [3, 13]) the solution of Eq. (7.2) is missed. Formally speaking, this 
solution originates from the singularities of the BFKL kernel. In practical terms, it means that 
we should be very careful approximating the full Balitsky-Kovchegov equation, by the simplified 
model of the BFKL kernel, namely, by diffusion approximation, in which this non-linear equation 
belongs to the Fisher-Kolmogorov-Petrovsky-Piscounov -type of equation (see Ref. [13]). This 
issue needs further careful investigation. 
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